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Mode-coupling theory (MCT) of sheared dense granular liquids is formulated. Starting 
from the Liouville equation of granular particles, the generalized Langevin equation is de- 
rived with the aid of the projection operator technique. The MCT equation for the density 
correlation function obtained from the generalized Langevin equation is almost equivalent 
to MCT equation for elastic particles under the shear. It is found that there should be the 
plateau in the density correlation function. 

§1. Introduction 

The rheology of granular materials is one of central concerns in granular physics. 
Although the granular particles exhibit unusual behaviors,^' Liu and NageP' sug- 
gested the possibility that the jamming transition of granular particles under the 
shear may be regarded as the glass transition at zero temperature. Since then the 
jamming has been recognized as one of the most important concepts in nonequi- 
librium rheology.'3®''3''3''3''Sl®'[iO},[i3,[ll,[l3} However, so far there exist few the- 
oretical works which discuss direct connections between the jamming transition of 
dense granular materials and those of the other materials such as colloidal suspen- 
sions^'^ from microscopic point of view. 

The lack of the direct evidences of the universality is due to the difficulties 
in describing the dynamics of dense granular media. On the other hand, indirect 
evidences of universality of the jamming transitions are being accumulated with 
time; (i) There are strong similarities for the emergence of dynamical yield stress in 
sheared granular material^^ to that of colloidal suspensions.'^ (ii) The steep peak 
of the nonlinear susceptibility so called X4 is observed in the vicinity of the jamming 
transition in experiments of shaken granular materials.'^ (iii) The long-tails in 
nearly elastic sheared granular fluids are almost the same as those for a system of 
elastic particles with a thermostat.'!^ 

We recognize that the gas kinetic theory can successfully describe the behaviors 
of relatively dilute granular flows including correlation effects,1^''2Sl'''SI>l22|i,l231i,l21},[23 
but it is obvious that we cannot use the kinetic theory in describing the jamming 
transition. On the other hand, many of theoretical work^''^''^''^ starting from 
statics of granular particles focus on the introduction of the effective temperature 
based on the idea of Edwards' compactivity in thermodynamics of granular media.f® 
Thus, so far there are few liquid theories in characterizing dense granular liquids. It 
should be noted that we do not have to introduce any exotic temperature such as 
compactivity in describing sheared granular liquids, because it is standard to use the 
granular temperature defined by T = — u)'^)/Nd where m, Vi, u, N and d 
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are respectively the mass of each particle, the velocity of i—th. particle, the velocity 
of the flow field, the number of particles and the dimension. In addition, it is known 
that the steady uniform shear flow has the uniform granular temperature except for 
the boundary layeif^'H'' in sheared dense granular systems, which is in contrast to 
the case of the relatively dilute sheared granular flows driven by the boundary.'^SJ 
Therefore, it may be possible to construct a liquid theory or the mode-coupling 
theory (MCT) by using the granular temperature to characterize the liquids near 
the jamming transition. 

In this paper we derive MCT equation for dense granular liquids from the Liou- 
ville equation of granular fluids. The derivation is formal and analogous to that for 
the ideal glass transition.l^S'ESliiES ^he organization of this paper is as follows. In 
the next section, we introduce the Liouville equation for granular gases. In section 
3, we formally derive the generalized Langevin equation by using the technique of 
the projection operator. In section 4, we apply the generalized Langevin equation 
for the field variables such as the density field under the shear. In section 5, we will 
derive the generalized Langevin equation for the density correlation function which 
is not a closed equation because of the formal representation of the memory ker- 
nel. In section 6, we introduce the mode-coupling approximation to obtain a closed 
equation for the density correlation function. The result is similar to that for ideal 
glass transition. In section 7, we discuss what we can predict from MCT and future 
problems. We also summarize our results. In Appendix A, we show the details of 
the derivation of the generalized Langevin equation. In Appendix B, we summarize 
the expression of the correlation functions in the presence of the shear. In Appendix 
C, we discuss the derivation of MCT for the sheared dense granular liquids. 



§2. Liouville equation 



Let us consider a system of identical hard spherical and smooth particles 
with their diameters a and the constant of the restitution e less than unity in 
d— dimensional space. The basic equation to describe the dynamics of the classical 
particles is the Liouville equation.l^'EB The Liouville equation for granular fluids 
has been discussed by Brey et q/.ES'EHI'EIIi some times ago. Let iCtot be the total 
Liouvillian which operates any physical function A{r(t)) as 

^^^£pi = i/:,,,A{r{t)), ^(r(t)) = e'^-*^(r(o)), (2-i) 

at 

where r{t) is the phase variable. The total Liouvillian consists of three parts, the 
free part, the collision part and the shear part. We denote iCtot as iCtot = iC + iCs, 
where iC is the sum of the free part and the collision part 

N 



iC 



E^^-T^ + ^EE^:'-^ (2-2) 



5r,- 2 



with the velocity Vj of j— th particle and the collision operator 

Tjk = o'''^ J d&5{rju - cT)0i-g,k ■ • 9jk\ib,k - 1), (2-3) 



Mode- coupling theory of sheared granular liquids 



3 



where rjk = rj — r^, gjj. = Vj — v^, 0{x) = 1 for x > and 0(x) = for otherwise. 
<T is the unit normal vector at contact and cr = a&. Here hjk is the cohision operator 
acting on any function X{vj^ViS) as 

bjkX{vj,Vk) = X{bjkVj,bjkVk) = X{v'j,v',,) (2-4) 

where the precolhsional velocities vj and change into the postcollisional velocities 
v'j and i)^, respectively. When bjk acts on gf^-j,, gjf^ changes as 

bjkQjk = 9jk - (1 + e)(a- • Qjk)^. (2-5) 
The Liouville operator to represent the shear flow iCg is 

^ ' 5 d 



'^^ = '^L.[y^9^^--y^^d^J^ (2-6) 

where the macroscopic velocity under the shear flow is assumed to be 

up = jyd/s^r,, (2-7) 

which recovers the equation of motion under the shear flow.'^ We also note that 
eq. (|2-6p can be derived from the DOLLS Hamiltonian.lSB'ES) 

It should be noted that the Liouvillian is not self-adjoint because of the violation 
of time reversal symmetry for each collision. The adjoint Liouvillian is defined for 
the equation of the phase function or A^— body distribution function p{r{t)) 

Pirn) = e-^-V(r(0)), = -iCtotpim). (2-8) 

The average of a physical quantity is defined as 

{Ait)) ^ J drpinAim) = J drA{r)p{r{t)). (2-9) 

From (j2Tp . (j2-8p and (j2-9p we obtain the following relations 

J drp{r)e'^""^A{r) = J drA{r)e-'^""^p{r) (2-10) 

and 

J drp{r)iCtotA{r) = - j drA{r)iCtotp{r). (2-11) 

The adjoint Liouvillian is obtained from the integration by parts of eq. (12-lip l^ 



N d 1 

iCtot = E • + 2 E E ^^^'^ + (2-12) 



where 



f,k = a''-^ J daO{& ■ gjj,)\& ■ g.j^\ [e-H{r,k - cr)b-^ + 6{rjk + ct)] (2-13) 
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and 

Here bjj^ is the inverse operator of bjk which satisfies bjj^bjk = bjkbj/^ = 1 and 

bjk9jk = 9jk - ^-J^i^ • 9jk)^- (2-15) 

Before closing this section, we add some remarks. First, since iCtot is not self- 
adjoint, the evloution operator acting on A*{r{t)) which is the complex conjugate 
of A{r{t)) is neither —iCtot nor —iCtot in (|2-8p . Second, the collision operators Tij 
and Tjj are reduced to the known Liouvillian for hard core particles in the limit of 
e = iW^ Third, the collision operators are exact for inelastic hard spherical granular 
particles for any density. 



§3. Generalized Langevin equation 



What we are interested in is the time correlation function between A{t) = 
A{r{t)) and the abbreviation of B{r{{})), as 

CAB{t) = {A{t)B*) 

= j dr[e'^'°*^A{r)]B*{r)p{r) = J drA{r)e-'^'°*^[p{r)B*{r)], (3-i) 

where we use the definition of iCtot for the last expression. The final expression of 
eq. (|3-l|) is obtained from the integration by parts with the assumption of zero net 
currents through the boundaries. 

It is not easy to handle eq. (j3Tp . because the adjoint Liouvillian iCtot has the 
propertji^ 

iCtoMr)B*{r)] = iCtot[p{r)]B*{r) + p{r)iC;,,B*{r), (3-2) 

where 

z£ro* - E ^ - 2 E E ^"fc + (3-3) 

j=l ^ j k=ij 

with 

Tjk ^ / da0{a ■ g^,)\a ■ g^kl^ir^k - <T){b]^ - !)• (3-4) 

However, when the initial phase function p{r) is invariant in the time evolution, i.e. 

Ctotpin = 0, (3-5) 

the simple treatment can be used. This situation may be realized in the case of 
the steady shear problem. In the later part of this paper, we only discuss the cases 
satisfying eq. ()3-5p . Note that {A{t)) is time- independent, though A{r{t)) can be 
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time-dependent under the steady condition (j3-5p . When the steady condition (j3-5p 
is satisfied, CAsit) is reduced to 

CAB{t) = I drp{r)A{r)e-'^^ott[B*{r)] = {AB*{-t)), (3-6) 

where we use the translational invariance of the correlation function, i.e., CAsit) = 
{A{t)B*) = {AB*{-t)). Thus, from eq.dMD with the aid of eq.(I33]) we obtain 



A*{r{t)) = e^^-U*(r); ^^^^£^ = ^CiotA*{^{t)) (3-7) 



by the replacement of B*{r{t)) by A*{r{t)). 

We may arise the naive question whether it is possible to derive the generalized 
Langevin equation for sheared granular fluids under the assumption of eq. (j3-5p . The 
answer of this question is " yes" . 

The procedure to derive the generalized Langevin equation is parallel to that of 
the classical simple liquid^ except for the non-self adjoint properties of iCtot- 
The details of the derivation are given in Appendix A. 

The final expression of the generalized Langevin equation is 

A{t)-iQA{t)+ [ dsM{t- s)A{s) = R{t). (3-8) 
Jo 

Here the memory kernel M{t) is given by 

M(*) . », (3.0, 



and i? is 



The random force in eq. (j3-8p is 

Rit) = e'^'^'^'^'QCtotA, R = iQCi^tA, (3-11) 
where we introduce the projection operator 

'^B{t)^ ^^^ff^ A; Q^l-V, (3-12) 

and the inner product 

iA,B{t)) = {B{t)A*) = j drB{t)A*p{r). (3-13) 

It is one of the most important results in this paper to obtain the generalized 
Langevin equation (I3-8P supplemented by eqs. (|3-9j) - (13-13p . Although nobody has 
discussed the generalized Langevin equation for granular fluids, the equation is useful 
in particular to construct a "liquid theory" of granular fluids. 
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§4. Some formulae in the presence of the shear 

The treatment of the generahzed Langevin equation (j3-8p for sheared granular 
liquids is still difficult because of the existence of iCs, but the effects of shear can 
be absorbed for the description of the hydrodynamic variables. For example, the 
operation of iCs to the local density field n(r, t) which is defined by 

N 

nir,t) = Y,^{r-rj{t)) (4-1) 

i=i 

for the position of j—th particle rj{t) is given by 

d d 



dxj ^ ^ dxi 

d ^ d 



~^y^^^^'^ ~ '^J^*)) = -iy-Q^ri{r,t). (4-2) 



3 

Therefore, the density field obeys 



+ iy^n{T, t) = iCn{r, t). (4-3) 
ot ox 

Similarly, any field variable A{r,t) obeys 

dA(r,t) dA(r,t) d , ^ ^ ^ 

-±^^^ + A(r.t)=.CA(r,t), (4-4) 

where iC is the sum of the free motion and the inelastic collisions acting on the 
hydrodynamic variable A{r,t). Thus, the time evolution in sheared systems is gov- 
erned by the Liouvillian iC instead of the total Liouvillian iCtot- Therefore, iCtot and 
iC^ot previous section and Appendix A can be formally replaced by iC and 

iC~ in the later discussion, respectively. We also note that our problem is reduced 
to a standard setup of field variables in the presence of the shear .13 'ISSI'lIS) 

It should be noted that the second term in the left hand side of eq. (|4-4p can be 
eliminated in the sheared frame 

A{rt, i) = A{r, t), n = r - ^yte^, (4-5) 

where t = t and is the unit vector along x— coordinate in the experimental frame. 
Indeed, A(rt,t) obeys 

dA{rt,i) 
di 

where we use 



iCA{rt,t) (4-6) 



dt = di - jydi, dy = dy- ^td^, d^^ = d^^ for a = 1, 3 (4-7) 



in the three dimensional case. 
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Let us introduce the Fourier transform 

Aq{t) = j dre'l-'-Air^t), Aq^{i) = J drte'^^''^ A{rt,i) (4-8) 
where is the stretched wave number related to the simple wave number q = 

(.qx,qy,qz) 

Qt^ {'lx,qy + iqxt,qz) (4-9) 

in the three dimensional system. It is easy to show the equivalences of these two 
representations 

Aq^{i)=Aq{t). (4-10) 

We also obtain the time evolution of the field variable Aq^{t) in the sheared frame 
is given by 

^^Aq^{i)=iCqAq^{i), (4.11) 

where we use 

dt=di + ^q,d,^, A = J_+^J_, A = ^ for a = 2, 3. (4-12) 

The right hand side of eq. (|4-lip includes the Liouvillian iCq which is the result of 
the Fourier transform in the sheared frame. 

Here, we note that the introduction of Fourier transform implicitly assumes to 
use the periodic boundary condition such as Lees-Edwards condition. Fortunately, 
it is known that rheological properties of dense granular flows under the shear are 
little affected by the choice of the boundary conditions .l^'l^ Thus, we believe that 
our theoretical argument in this paper can be used even in physical situations. 

In Appendix B, we summarize the relevant expression of the correlation function 
in the presence of the shear. This Appendix might be useful to resolve the confusion 
among various expressions in literature. 

§5. Generalized Langevin equation for the sheared granular liquids 

In section 3, we have derived the generalized Langevin equation, but the result 
is too formal to apply it to physical processes. Since our objective is to derive the 
mode-coupling theory (MCT) for sheared granular liquids, we had better focus our 
attention on the density fluctuations. For this purpose we first summarize what we 
calculate to derive MCT equation. 

MCT is the theory to describe the time correlation function of the density fluc- 
tuations. The density of liquids is given by eq. (|4-ip and its Fourier component in 
the experimental frame becomes 

N 

nq{t) = Y^e^^-^^^'\ (5-1) 

i=i 

From the argument in Appendix B, the density correlation function in the experi- 
mental frame should be characterized by F(q,t) defined in eq. ()B-5P where suffices 
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A = B for the density field in Appendix B are eliminated. Let us use t) for 
later discussion, in which the explicit definition is given by 

F{q,t) ^ l{nqjt)n^q{0)) = 1 ^exp[i(q_, • rj{t) - q ■ r^)]. (5-2) 

The function F{q,t) is reduced to the scattering function at t = 

F{q,t = 0) = l(ng(0)n_q(0)) ^ S{q). (5-3) 

Note that the form of the structure factor S{q) or the pair correlation function g{r) 
which is the Fourier transform of S{q) for granular fluids has not been established. 
There are some theoretical studies^'^' without any external forces but no theoretical 
work in the presence of the shear. However, the theoretical approaches suggest that 
the structure factor is independent of the relaxation processes in granular liquids, 
which should be valid in our steady case ()3-5p . Therefore, we can separate the 
problem to determine F{q,t) from the determination of S{q). We also note that 
the anisotropy of the pair correlation function induced by the shear exists but is 
small in the simulations of sheared granular flow^^'HS as in the case of sheared 
colloidal particles. One of the characteristics, however, for granular liquids appears 
in g{r) where the first peak around \r\ = o" is higher than that in the conventional 
cases.HB'HS 

Consider the density defined in eq. (j5-l|) . Its time evolution is described by 



d .... ,L 



dt 



nqjt)=iq.tjijjt). (5-4) 



Here we introduce jg ^ (t) = q ^ ■ jg ^ (t) with the current 

h-M^^Y.^^i^y'^-'"''^'^ (5-5) 



and q^t = q_t/q-t- 

To describe the slow process of the density correlation function in sheared gran- 
ular liquids we formally replace A{t) and iCtot in section 3 by 




and i^q_^, respectively, where 5ng_^{t) = e*^-*''"^^*-' — (27r)'^n(5(q_j) with the 
average density n, and the longitudinal mode of Fourier component of the current 
field. 

Therefore, eq. (|3-8p is replaced by 



dAgJt) 
dt 



iOg ■ Agjt) - ^ dsMg_^^_^^{t - s) • Ag_^{s) + Rgjt), 



(5-7) 
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where the memory kernel, the random force and the characteristic frequency are 
respectively given by 

Mqjt) = {Rq,RgJt)){Ag,Aq)-\ (5-8) 

Rq (t)=eMi f QCq Mdt']iQCqAq, Rq=iQ£qAq, (5-9) 
Jo 

and 

if2q = {Aq, Aq){Aq,Aq)-\ (5-10) 
We also introduce the correlation matrix 

Cq jt) ^ {AqJt)A*qm = {Aq,Aq Jt)). (5-11) 
Using {A, R{t)) = we obtain 
dCq_^{t) 
dt 

where Cq_^{t) is 



iOq ■ Cq_^{t) - ^ dsMq^^^_^^{t - s) ■ Cq_^{s), (5-12) 



( {6n.qr^q_^{t)) {6n.qj^Jt)) \ 

which reduces to 

C,(0) = ( "^„<«> Sr) (5.14) 



at the equal time case (t = 0), where we use the granular temperature as T. Here, 
we should add some explanations for eq. (|5-14p . The diagonal elements of (j5-14p are 
respectively the definitions of the structure factor and the granular temperature. 
On the other hand, the off-diagonal element becomes {j-qSnq) = "^Zii^ ' Pi) ~ 
m^^ ^•((q-Pj)e~*^'^'(5(g))(27r)^n, where p,j = m{vi—^yex) is the linear momentum 
of the particle i. It is natural to assume that the momentum is independent of the 
position and the average of linear term of the momentum is zero in the statistical 
average of the uniform shear. Thus, we assume that the off-diagonal element in 
eq. (|5-14p is zero. 

From eqs. (|5T0p . (j5-lip . (j5-13p and (|5-14p we also have the expression 

Here we assume that the correlation between a field variable and its time derivative 
is always zero. We also use the following identity 

= -^yZ'y^i'^j-^y^^f) = ^—^ (5-16) 
md ^ m 
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to obtain eq. (|5-15p . It should be noted that the granular temperature T is propor- 
tional to 7^ in granular fluids under the steady shear.^^ For simplicity, we assume 
that the temperature is uniform, which can be realized in a periodic system in the 
vicinity of jamming transitions. The uniform temperature is also realized in the bulk 
region of dense granular flows under the physical boundary. 

The complex conjugate of the random force at t = is given by 

R*q{0) = iQ*CZqA*q = i{l - m^Zq ( ^Jl^ ) , (5-17) 

where V* and Q* are the complex conjugate of the projection operators V and Q, 
respectively. From eq. (|3-7p with the replacement of by £g we obtain the relation 



The second term in the last expression of eq. (|5-17p can be rewritten as 

^^^g^g- ^A*q,A*q) 

Therefore, we obtain 



Ra=\ dj-q , I = ( ). (5.20) 



Similarly, Rq_^ (t) is given by 

where we use eq. ()5-9p . 

Let us look at the equation of motion term by term. First, the left hand side of 
eq. (|5-12p becomes 



dCqJt) _ ( f^{6n.q6nqjt)) i{6n^qj{jjt)) 



dt I iij'^q^nqjt) Hj'lqjijjt)) 



(5-22) 



in our problem by using eq. (j5-13p . Note that the lower left corner element becomes 
{N/iq){d?F{q,t)/dt^). Second, from eqs. ([FT3]) and ([FT5]) we obtain 

,„ r (f\ ( ° n\( i^n^q^nqjt)) {6n^qJ^Jt)) \ 
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where the lower left corner term is —{qNT/imS{q))F{q,t). Lastly, the memory 
matrix is 

/ n n \ ( ^ \ / 

^1-M=[ iR_qRq ^(t)) ) ""f^ ^ = "^^^-^^^-(^^^ 

\ NT / ^ NT 

(5-24) 

where we use eqs. (|5-15p and (j5-20p . Concerning the lower left corner, using eq. (j5-12p 
we obtain 

^ + Mr/''-'^ + ^ - = 0. (5.25) 

§6. Mode-coupling approximation 

The derivation of eq. ()5-25p is almost exact except for the assumption of the 
uniform temperature but eq. (l5-25p is not a closed equation of F{q,t). On the other 
hand, we are only interested in the slow dynamics of F{q, t) in the vicinity of the jam- 
ming transition. One of the possible approaches to get a closed equation for -F(q, t) 
is to diagonalize the linearized hydrodynamic equations of Snq(t) and jq{t)^^^^ 
The quantitative validity of this approach for dilute granular gases has been con- 
firmed, but we may not use this approach in the dense granular liquids. As another 
approach we adopt the mode-coupling approximation for granular liquids near jam- 
ming transition in which the argument borrows from that for the conventional glass 
transition. We should stress that the derivation of MCT equation is not the final 
goal to describe the jamming transition. Indeed, we will need to determine S{q) and 
the quantitative relation between T and 7 for granular fluids. 

Let us consider the term {R^qRq (t)). This can be rewritten as (^_ge*^^^-*^*i?q) 
from eq. (|5-8p . This term may consist of the fast part and the slow part. The fast 
part may be approximated by rq_^^/T5{t) with a friction constant -Tq.t, since the fast 
friction may be proportional to the kinetic viscosity and in the hydrodynamic 
limit. ^^-^ The slow part represents the important contribution for the slow relaxation 
process, which is nothing but the mode-coupling memory kernel. 

(i) Let us replace the slow part of e by V2, where we introduce 
the new projection operator acting on any function B 

-p^B^ ^Q....9.-.(^g3-.Q.-.^)^^g.-.9../g3-.9.-.)"^ (6-1) 

Ql tCl2 vQs tQ4 t 

where Aq g = 5nq Snq . The projection operator V2 is the simple projec- 
tion onto its dominant slow mode. We neglect the contribution from Q as in the 
case of the conventional MCT. 

(ii) We assume the factorization of four-point correlation into the product of 
two-point ones to obtain a closed equation for the density correlation function. 

Using these approximations the derivation of MCT equation is straightforward. 
The details of the derivation are summarized in Appendix C. 
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Finally, we summarize MCT equation 

(6.2) 

with 

<y^ = T6k^l ^^-*l^q-fc,fcV.-fc-./^-J^(^^*'*)^(^-*-^-*'*) (6-3) 

for d = 3, where are connected with q by eq. (j4-9p . where Vg ^ is given by ()C-14p 
with the aid of the direct correlation function (jC-lOp . As suggested in section 5, the 
difference between sheared dense granular liquids and the conventional cases of the 
sheared MCT equations appears through (i) the forms of the structure factor and 
the direct correlation function, and (ii) the granular temperature which disappears 
in the limit of the low shear rate. Although the second term in the left hand side in 
(j6-2p looks dominant in the low temperature limit or the low shear rate limit at the 
first glance, all the terms are in the same order when we assume Bagnold's scaling^SJ 
T ~ 7^ and the time is scaled by 7""'^. 



§7. Discussion and conclusion 



7.1. What can we predict from MCT equation for sheared granular liquids? 

In this paper, we have demonstrated that the MCT equation can be derived in 
sheared dense granular liquids. This is the first step to understand the universal 
feature of the dense granular liquids and the jamming transition. Let us summarize 
what we can predict for sheared dense granular liquids from MCT equation (j6-2p 
supplemented with (j6-3p . 

First, we expect that there is a plateau in the relaxation of F{q,t) at least for 
nearly elastic granular particles when the density is larger than a certain threshold 
value. Indeed, eqs. (j6-2p and (j6-3p reduce to the well-known MCT equation for rela- 
tively short time t ^ r(7, e) where the characteristic time (or the life time) of the 
plateau t(7, e) may satisfy r(7, e) oc 7""*^ at least for nearly elastic cases.ll^ Then 
MCT equation exhibits a quasi-arrested state of particles in cages. However, this 
arrested state is destructed by the stretching of cages induced by the sheared force. 
Therefore, we may write 

F{q, t) ^ ^(T(7, e))F{q, t/r(7, e)) (7-1) 

in the quasi-arrested state. Here <P{T) represents the scale factor as a function of the 
temperature, which tends to unity in the limit of T ^ or 7 ^ 0. We also note that 
there is no steady state in the limit of e = 1 because of the viscous heating effect .33 
Although so far there was no report of the existence of visible plateau in granular 
liquids,"^ we have reproduced a two-step relaxation of F{q, t) from the simulation 
of a dense and nearly elastic sheared granular liquid as explained later. Similar to 
the conventional cases^- to reproduce the two-step relaxation, we need to prepare 
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binary systems. Indeed, there is the crystalization for mono-disperse spheres, while 
there is no plateau for randomly dispersed particles. In addition, we note that the 
range of parameters to observe two-step relaxations seems to be narrower than the 
conventional cases. 

Let us explain the preliminary result of our simulation briefly. The system 
simulated is a three-dimensional 80:20 mixture of AT = 1000 Lennard-Jones system 
in which the potential is given by 



a/3 



Tap 



12 



Tap 



(7-2) 



where a and (3 refer to two speicies, and raf^ denotes the distance between the 
particle a and the particle (5. The particles arc confined in a periodic box whose 
linear dimension is dAaAA under the Lees-Edwards boundary condition. We choose 
GAB = ^-^^AA, ^BB = 0.5eAA, asB = 0.88cJAyi and eab = 0.8c7aa- We introduce the 
dissipative force -r]{vaf3 ■ r af3)r / r^f^ with r^is = ra-rp and Va0 = Va-vp for 
r^p < 1.12246cj„^ in the equation of motion of the particle a, where the dissipation 
parameter rj = O.OOl in the dimensionless unit. Note that all quantities are non- 
dimensionalized by the particles diameter aAA, the interaction energy caa and the 
time To = {mAcr\j^^/ caa)^^"^ with the mass of A particle tha- We add the shear with 
7 = O.OOl/ro to the system. Figure 1 is the result of the mean sqaure displacement 
of particles {r{t)'^) = — ^i(O))^)/-^) where Zi{t) is the position of i— th 

particle in the vertical direction to the sheared plane. It is clear that the particles 
are in a quasi-arrested state in the middle stage of time. Figure 2 is the result of 
F'^{q,t) = X]j(cos(g(zj(t) — Zi{0))) /N with g = 15 in the dimensionless unit. We find 
that F{q, t) has the plateau in the middle of the relaxation process. The details of 
our simulation will be reported elsewhere. 



<r(t)^ 




0.01 



0.0001 



0.01 0.1 1 10 100 1000 

t 

Fig. 1. The mean square displacment of particles {r{t)^) = X^j{(-?<(t) — Zi{Q)Y) /N as a function of 
time . 



Second, it is difficult to describe the jamming transition based on the hard-core 
model which we use in this paper. Let us demonstrate this difficulty as follows. 
We expect that the plateau r(7, e) oc becomes long as 7 decreases. Eventually, 
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0.01 0.1 1 10 100 1000 
t 



Fig. 2. The relaxation of the correlation function F"{q,t) — '^i{cos{q{zi{t) — Zi{0)))/N. 



the plateau becomes infinitely long in the limit of 7 ^ 0. This is the reflection 
of no motion of particles in the limit of 7 0. When we adopt the factorization 
approximation as in the framework of MCT, the shear stress axy for sheared granular 
liquids can be represented by a function of S{q), S{qf) and F{q,t) as^^' 



where S'{k) = dS{k)/dk, and we ignore anisotropy of the structure factor and the 
density correlation function. To derive eq. (|7-3p we use the decoupling approximation 
and the Green-Kubo formula which may be suspicious in granular liquids. However, 
the decoupling approximation is consistent with MCT and the Green-Kubo formula 
gives, at least, a good approximate expression for the nearly elastic granular liquids. 
Thus, the expression (l7-3p should be valid for the nearly elastic granular liquids. If 
we assume that there is no relaxation of F{k,t) in a quasi-arrested state, eq. (|7-3p 
may be replaced by 

,. ,,2 r^^'^^ , r „ ,4S'{k)s'ik) 

oc T^F{k, r(7, e))V(7, e) oc T, (7-4) 

where we replace the integrand by that for t <^ r(7,e) ~ 7"^- From eq. (j7-4p the 
viscosity r] can be evaluated as 77 ~ T'ItI"^ as in the case of colloidal suspensions 
near the glass transition.lIS'lSli'lISJ'llB This does not mean that the sheared granular 
liquid displays the shear thinning property, but this result is consistent with Bag- 
nold's scaling ~ I7I or T ~ 7^ for hard-core granular fluidd^ which is one of the 
shear thickenning relations. Indeed, if we assume that the continuity equation for 
the energy is still relevant, the uniform steady granular liquids obey rfj'^ ~ T^/"^, 
Therefore, we obtain the relation 

T ~ 7^ (7-5) 



which is nothing but Bagnold's scaling. 
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From eqs. (l7-4p and ()7-5l) we may obtain 

<7xj/ ~ T ~ 7^. (7-6) 

We should note that Bagnold's scahng ()7-5p is violated in the vicinity of the jamming 
transition as T ~ 7^ with /? ~ 1.1. SS) -^ye qI^q note that the jamming may be the 
continuous transition 39 '33 which is different from the sheared dynamic yield stress 
at a constant temperature. 

Equations (j7-6p may imply the following two things: (i) The stress can be in- 
depent of the shear rate when we keep a constant temperature with changing the 
restitution constant e. (ii) However, the stress becomes zero in the limit of 7 — > 0. 
Thus, we may conclude that our MCT for hard-spheres is invalid to describe the jam- 
ming transition. This conclusion is reasonable, because the jamming transition is 
originally defined by the point of non-zero bulk and shear moduli at zero temperature 
without kinetic energy.SSJ 

7.2. Future problems 

In this paper, we only present the formal derivation of MCT equation and give 
some qualitative predictions. To know the quantitative details of the behaviors of 
dense granular liquids, we need to know the explicit expression of S{q) and the 
relation between the granular temperature and the shear rate. As long as we know, 
we do not have any satisfactory theory to describe sheared dense granular liquids. 
This will be our important subject to the future. 

We should note that our MCT is based on the mean field description in which 
the system is almost uniform. The heterogeneity may be important in the actual 
jamming transitions, but this heterogeneity may disappear if we only observe two- 
point correlation functions as in the case of the glass transition. We believe that 
we will need to formulate the model including higher correlation to describe such 
heterogeneity. 

We also note that our starting equation is the Liouville equation for hard spher- 
ical particles. The system of inelastic hard spheres may cause the inelastic collapse 
in which the collision frequency is divergent when the density becomes extremely 
high in the vicinity of the jamming transition.'^' We also note that the jamming is 
caused by the simultaneous contacts among particle^^ which cannot be described 
by hard-core models. Although we obtain the formal expression for MCT equation, 
we may have to extend our formulation to soft spheres to describe the jammed state, 
which was not obtained previously. 

In this paper, we adopt the formal derivation based on the projection operator 
formalism. The advantage of this method is to get the exact expression in the mid- 
dle of calculation, in which equation (j5-25p should be exact. On the other hand, the 
procedure to obtain MCT equation in section 6 contains some uncontrolled approx- 
imations. Therefore, it is difficult to systematically improve the decoupling approx- 
imation used to derive MCT equation. There are several efforts to overcome such 
difficulties in conventional glass transitions^ ''^''^''^ based on the field theoretic 
approaches. Such methods may need to improve our treatment presented here. 
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7.3. Conclusion 

We summarize what we have carried out in this paper, (i) We have derived the 
generahzed Langevin equation (l3-8p with (l3-9p for granular fluids under the steady 
condition (j3-5p . This equation is formal and the derivation is exact, (ii) We apply 
the generahzed Langevin equation to the sheared case in which the stretching of the 
wave number in Liouvihian is enough to describe the sheared system in eq. (|4Tip . 
This equation is also believed to be exact, (iii) We derive the equation for the density 
correlation function F{q,t) as in eq. (|5-25p . The equation is almost exact under the 
assumption that the granular temperature is uniform, (iv) We derive MCT equation 
(]6-2p with (j6-3p in which we include some uncontrolled approximations as we have 
used to derive the MCT equation for conventional glassy systems, (v) We suggest 
the existence of a plateau of F{q, t) for the dense granular liquids, (vi) We also 
indicate that hard-core models are insufficient to describe the jamming transition. 
It should be noted that the MCT equation includes the static structure factor of 
sheared granular materials which should be different from the conventional cases. 
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Appendix A 

The derivation of the generalized Langevin equation 13- 8\) 

In this Appendix we demonstrate how to derive the generalized Langevin equa- 
tion (l3-8p . The procedure is parallel to that for the simple liquids.'^''^ 
Let us introduce 

With the aid of the Laplace transform 

/•OD 

Aiz) = / dte"^A{t), (A-2) 



eq. (|2-ip is reduced to 

A{z) = {z + Ctot)-^iA{0) (A-3) 
Operating V and Q to eq. ()A-3P we obtain 

zrA{z) + VCtotVA{z) + VCtotQA{z) = iA (A-4) 
zQA{z) + QCtotVA{z) + QCtotQA{z) = 0. (A-5) 
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Substituting eq. (IA-5p or the equivalent form QA{z) = —{z-\- QCtotQ)^^ QCtotP ^{z) 
into eq. (|A-4p we obtain 

zVA{z) + VCtotVA{z) - VCtotiz + QCtotQT^ QCtotV A{z) = iA. (A-6) 

Using the relation {A,VB) = (A, B) the inner product of eq. (|A-6p with A becomes 

z{A, A{z)) + {A, CtotVA{z)) - {A, Ctot{z + QCtotQ)-^ QCtotV A{z)) = i{A, A). (A-7) 

From the relation 

{A.iCtotB) = {{iCtotB)A*) = J dr{iCtotB)A*p{r) = - j dr{iCt,^A*)Bp{r) 

= -{{zC^,,A*)B) = -{iC^,,A,B) (A-8) 

eq. ()A-7p is reduced to 

- {A,CtotVA{z)) {Ct,^A,{z + QCtotQ)-^QCtotVA{z)) _^ 

{A, A) + {A, A) ^ 

where Y{z) is the Laplace transform of y{t). 

Here the second term in the left hand side of eq. ()A-9P can be rewritten as QY{z) 
with A = dA/dt at t = 0. To derive 17 in eq. (j3-10l) we use the relation 

{A,CtotVA{z)) _ {A^CtotA) {A,A{z)) _1{A,A)~^^ , 

{A, A) ~ {A, A) {A, A) ~i{A,A) ^ '- ^ ' 

On the other hand, the numerator of the third term in the left hand side of ()A-9P 
can be written as 

(A-ot^, {z + QCtotQ)-'QCtotVA{z)) = {Ci^^A, Q{z + QCtotQr'QCtotA)^-^^^^ 

= {C^^^A, Q{z + QCtotQ)-^QCtotA)Y{z) 
= {QCi,^A, {z + QCtotQ)-^QCtotA)Y{z) 
= -{R,{z + QCtotQ)-'R)Yiz), (A-11) 

where R is defined in eq. (|3-lip and 

R = iQCtotA., R = iQC-t^tA. (A-12) 

Here we use the relations 

{QB, C) = ((1 - V)B, C) = {B, C) - |^(^, C) (A.13) 

and 

(B, QC) = {B, (1 - V)C) = {B, C) - C) (A-M) 

with {A, B) = (B, A) for any real functions A, B and C. (For complex functions, 
the relation {B,A) = {A*,B*) should be considered). 
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Introducing the memory kernel in the Laplace form 

M(z) = i{R, (z + QCtotQy^R)iA, Ay\ (A-15) 
eq. ()A-9p can be rewritten as 

-i{z + n)Y{z)+M{z)Y{z) = l. (A-16) 

Therefore, we obtain 

Y{t)-iQY{t)+ [ dsM{t-s)Y{s) = (A-17) 
Jo 

in terms of the inverse Laplace transform, where M(t) is the memory kernel. 

Equation ()A-17p describes the time evolution of the correlation function. On the 
other hand, the time evolution of A{t) is affected by the fluctuating force. Indeed 

A^{z) = QA{z) (A-18) 

obeys 

{z + QCtotQ)A'^{z) = -QLtotVA^iz) = -QCtotY{z)A 

= --RY{z), (A-19) 
i 

where we use the definition of Y{z), eqs. ()A-5P and ()A-12p . Thus, we obtain 

A^{z) = Y{z)R{z), (A-20) 
where R{z) = i{z + QCtotQT^R- Substituting eo. dXTBl) into ea. dX^ we obtain 

{-iz-iQ + M{z))A^{z) = R{z). (A-21) 

Therefore we obtain 

A{t) - iQA{t) + / dsM{t - s)A{s) = R{t) (A-22) 
Jo 

with 

A{t) = QA{t), R{t) = exj>[iQ£totQt]R. (A-23) 
Since A{t) satisfies A{t) = YA{t) + A{t) we obtain 



Appendix B 

The expression of the correlation function in the presence of the shear 

Let us summarize the correlation function in the presence of the shear. The 
correlation function should satisfy the translational invariance condition 



{Airt + a,i)Bir' + a,0)) = t)5(r', 0)) (B-1) 
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where a is the shift vector acting on all particles in the sheared frame. The Fourier 
transform of both sides of this equation leads to 

CAq^Bg,{t) ^ {Aq^{t)Bq>m = e'^^^^^^l'^^ {Aq^{t)Bq,m . (6-2) 

Thus, we obtain 

CAq^BqXt) = {Aq^{i)Bq,{0))5q^^_q, (B-3) 
in the sheared frame. From the relation ()4-10p we can write 

CAq^B_q,{t) = {Aq{t)B_q,{0)) = 5q^,q>FAB{q,t) (B-4) 

where FAsiQjt) = {Aq{t)B-q^{0)) in the experimental frame. This relation is ob- 
tained by Fuchs and Cates.ll^ This relation can be rewritten as 

CAq,,B^q{t) = Sq',q_FAB{q,t) (B-5) 

where F^B(q,t) = pABiQ-i.t) = {Aq_^{t)B_q{0)) = {Aq^{t)B*q{0)) , which is ob- 
tained by Miyazaki et al!^ 

Appendix C 

The derivation of MCT equation 



Let us explain the details of the derivation of MCT equation in this Appendix. 
From the two approximations mentioned in section 6 we get 

V2RqJt)= Yl yq-Mi-v'l2-t)^nq,Ji)^nq^Jt), (C-1) 

where 

, ^- {5nq^Jt)5nq,^Jt)RqJt)) 

Vq_M.-.q,^t) - 2. l^8nqat)5r^qat)5nqat)5r^qat)y ^^"'^ 

The product of the four density fields in the denominator of eq. ()C-2p may be factor- 
ized into the products of two structure factors. The numerator of eq. (|C-2p is 



( Sn^Jt)5nf^_^_gJt)RqJt)) 

djq_ (t) 

. q-tT 



mS{q- 



{6n_^_^ {t)6nj^_^_q_^ {t)6nq^^ (t)). (C-3) 



The first term of the above equation is 

djq d d 

i^^kj^k.^^q ,-^) = -ij^^^-kj^k.,-qjqj-i^-kjj^^nk.,~qjjqj 

(C-4) 
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Let us calculate the first term of the right hand side of eq. ()C-4p . 

j k 



I 

T 
-i — 
m 



= -t{k^^.q_^)LNS{k^t-q-t), (C-5) 
m 

where we use {vj^a^i^p) = 5ji5apT /m. The other term similarly gives 

- i^^^kjt)jti^^k^,-qjt))jqjt)) = iiq-t ■ {q-t - k^t))^NS{k^t). (C-6) 

The last term of eq. ()C-3p is hard to compute directly. When we adopt trhe decoupling 
approximation or Kirkwood approximation, we obtairPBS 

{5n_f^ Jt)Snf^_^_g Jt)dnq Jt)) ^ NS{k^t)Siq_t)S{k^t - q-t)- (C-7) 

Since q„( satisfies the periodic boundary condition, we expect that the vertex func- 
tion Vq {ki_i,k2-t) is the function of the difference of the wave vectors k t = 
ki_i — k2-f Thus, we can write 

Vg_^(fcl-*,fc2-*) = ^fc„,,g_^„fc_,. (0-8) 

Subsitituting ()C-4p - ()C-7p into ()C-2p with taking the summation over k-f, we obtain 

T/ f q_^-k_t , q_t-{q-t-k-t) I . A 

uq^-k-, - ^ \ S{k-t) + Sik-t - q_,) ■ "^-''j 

%TtT 

= 1^ {{q-t ■ k-t)c{k_t) + q-t ■ {q-t - k-t)c{k_t - q-t)} (C-9) 
where we introduce the direct correlation function c(fcj^ as 

c{k) = -{l- S{k)-^). (C-10) 
n 

Similarly, we obtain 

V2R-q = Y.V^ q_k^n_j,5nt, (C-11) 

k 

where 

Vk,q-k = ^{(^ • k)c{k) + q {q- k)c{k - q)}. (0-12) 
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Therefore, we may obtain 



where 



k 



2m^ 



k-t,t) 



k-t,t), (C 



Vfc k = {i^- k)c{k) +q {q- k)c{q - k)}. 



(C-14) 



Using the equations obtained in this Appendix we obtain the final expression of 
MOT equation (HH) with (|6^ . 
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